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PAPERS 


TURBULENT TRANSFER MECHANISM AND SUS- 
PENDED SEDIMENT CLOSED CHANNELS 


Tests rectangular channel, 10.5 in. in. cross section, transporting 
water mixed with different proportions sand, are described this paper. 
terms mean sedimentation diameters, two sizes sand were used—0.10 
and The experimental data were used study (a) the effects the 
presence sand suspension the characteristics the flow, and (b) the 
momentum transfer coefficient and the concentration transfer coefficient the 
profile through the center line the channel, assuming two-dimensional flow. 

The present experiments show that the universal constant turbulent 
exchange, decreases with the increase the suspended material. The value 
decreases 0.20 when the average concentration 4.3% weight. 
limit was approached for the decrease The friction coefficient was hardly 
affected the presence the sand when the velocity flow was higher than 
certain rate. This velocity related the bed load; the velocity which 
all the dunes the bed are carried suspension; and increases with the total 
load the flume and the size the sand used. Below that velocity, exceeds 
that clear water. 

The sediment transfer coefficient was found equal 1.5 times the 
momentum transfer coefficient, for the 0.10-mm sand and 1.3 for the 
0.16-mm sand. Both and follow the normal parabolic form the outer 
two thirds the channel and have constant value the middle third. 


INTRODUCTION 
letter symbols introduced this paper are defined where 
they first appear, the text illustration, and are assembled alphabetically 
the Appendix for convenience reference. 
The modern developments the turbulent flow theory have created interest 
the sediment suspension and the heat flow problems. was Morrough 


comments are invited for publication; the last discussion should submitted 
August 1951. 


Engr., Ministry Public Works, Cairo, Egypt. 
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ASCE, who first made use the Taylor-Schmidt findings de- 
riving the equation for the distribution suspended sand turbulent flow. 


equating the rate upward transfer, suspended sediment 


resulting from turbulent exchange and the rate settling, under gravi- 
tational force, Mr. O’Brien obtained the equation: 


which measured vertically upward; the concentration suspended 
material; the settling velocity the sediment the still fluid; and the 
transfer coefficient for the sediment. 

Mr. von Hon. ASCE, gave the analogy between the transfer 
mass heat and the transfer momentum. showed that the co- 
efficient all these kinds transfer should have the same form, but stated 
that was subject further discussion whether those coefficients have 
the same values. 

Another step was taken Sherwood and who found 
that the transfer coefficient for water vapor turbulent gas stream was not 
equal the momentum transfer coefficient €m, but that they seemed bear 
constant relation each other. Other investigators disclosed that this 
approximation can assumed also case heat transfer. Therefore, may 
assumed that can related linearly 


which the coefficient proportionality. known that, for turbulent 
velocity; and the density the fluid. For two-dimensional flow, 
vertical distance the maximum velocity—that the plane zero shear. 


du/dy 


which the maximum velocity and Therefore, when 
du/dy evaluated from Eq. the result 


the Theory Turbulent Flow and Its Relation Mor- 
rough O’Brien, Transactions, Am. Geophysical Union, Vol. 14, 1933, pp. 487-491. 
Aspects the Turbulence Problem,” Theodor von Mechanical Engineering, Vol. 
57, 1935, pp. 407-412. 
Chemistry, Vol. 31, 1034. 
Mechanics Turbulent Boris Bakhmeteff, Princeton Univ. Press, Princeton, J., 


TRANSFER MECHANISM 


and, its dimensionless form, 


The universal constant used the left-hand side that the right-hand side 
will function only the position. 

Eqs. give zero values for both and the center the channel. 
Many investigators maintain that the transfer coefficients are not zero the 
center and, therefore, that the equations are not valid this point. 
explained how these relations were derived from the assumption 
similar flow pattern; then stated that the center part pipe the similar- 
ity assumption cannot correct. Brooks and Berggren’ showed how 
Prandtl tried correct that region and how the experiments Nikuradse 
gave definite values the center the pipes his experiments. They 
showed that the two reasonable assumptions which can made the center 
are either take constant according the results Messrs. Sherwood 
and Woertz‘ assume error curve for €m. 

constant over certain region, Eq. can integrated for that region 
give 


which the concentration any arbitrary reference level 
follows Eq. another region, Eq. will integrated give 


Vito ASCE, verified Eq. and showed how well fitted 
the experimental points open channel. also showed that the exist- 
ence suspended load tends suppress damp out the turbulence, causing 
decrease the value the von universal constant turbulent 
exchange The main lines Mr. Vanoni’s work were extended the 
present paper cover the study closed channel. 


and Skin Theodor von Journal the Aeronautical Sciences, Vol. 
1-2, 1934, pp. 1-20. 

Turbulent Transfer Across Planes Zero Momentum-Exchange,” Brooks and 
Berggren, Transactions, Am. Geophysical Union, Pt. VI, 1944, pp. 889-896. 
Suspended Sediment Water,” Vito Vanoni, Transactions, ASCE, Vol. 111, 

» Pp. 
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OBJECTIVE EXPERIMENTS 
The two main objectives the present experiments are follows: 


The first objective was continue the study the effects the sus- 
pended sediment the characteristics the flow. closed channel was 
selected for these studies instead open channel because the relative ease 
with which steady uniform conditions can established and accurate mea- 
surements made. attempt was made solve the partial differential equa- 
tion the second degree that represents the sediment distribution the case 
circular pipe, but the equation was too complicated solved reasonable 
time. the flow between two wide parallel plates was used. order 
study the effect the presence sand the wall roughness all the walls 
the flume were painted make them smooth possible. 

The second objective was study both the momentum transfer coeffi- 
cient and the sediment transfer coefficient their distribution over the 
cross section, especially near the center line the channel, and their interrela- 
tion. 

APPARATUS AND PROCEDURE 

Main experiments were conducted closed-circuit flume 
(Fig. about 10.5 in. wide, in. deep, and long. Precise measurements 
show that the section 2.986 in. 0.5%) 10.507 in. 0.4%). For the 
present experiments the flume was kept horizontal. 

The top cover the flume was made steel plate, supported from 
above, shown Fig. The cover was sealed the sides solid 
rubber rod in. diameter). All the surfaces were covered carefully with 
two coats bitumastic paint which gave reasonably smooth surface. 


(a) PLAN 


Direction Flow 


Transition 


Vari-drive 


Circulating Pump 


ELEVATION (c) TYPICAL 


The pitot tube and the samplers can traverse the cross section through 
opening the flume cover (Fig. 2). This working section 28.5 from the 
entrance the channel. the downstream side the working section, there 
were glass windows long each side wall and plastic window the cover 
10.5 in. long and in. wide. The water the mixture sand and water was 
circulated the flume propeller pump driven variable speed motor. 

Measuring standard Prandtl pitot-static tube in. dia- 
meter was used measure the velocity. The differential pressure the tube 
was read 0.001 water-air manometer and converted velocity 


=i is$Yf ! 
Motor 
Adjustable Venturi Meter 
Support 10.507 In. 
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assuming coefficient unity for the tube. Two samplers were used deter- 
mine the suspended sand concentration different points. One them 
reached the upper surface and the other reached the bed (see Fig. 
were made brass tubing with external diameter Their tips were 
flattened that the inside dimensions were 0.22 in. 
Each sample was siphoned from the flow 
and collected bottle 1-liter volume. 

The elevation the bottle and, hence, the head 
the siphon could adjusted give the same velocity 
the inlet the sampler the local stream velocity 
the sampling point. The effect the rate 
sampling the proportion sediment removed with 
the sample, well the sampling technique, has been 
clearly explained The sediment was re- 
moved from the sample filtering and was dried and 
weighed the nearest 0.005 

general sample was collected special sampler 
the end the flume. Three liters were collected for 


each run and they were used find average value 
concentration suspended sand over the section. 
This average concentration was called 


The pitot tube, one the two samplers, SECTION 


fixed carriage contrived that could move 
across and along the stream and the vertical di- 
rection. The three coordinates the instrument 


operator either locate the position the instrument 

The flow rate and the average velocity, 
channel were determined from velocity measurements 
with the pitot tube. From complete velocity trav- (b) PLAN 


the carriage could read from scales that enabled the 


erses runs and 71, the relation— 


—was found, which the average velocity 
along the center-line profile. all runs, except the 
aforementioned, the average velocity was determined 
this empirical relation. This method was 


adopted preference that using the venturi 
meter. The meter behaved erratically, probably 


because the presence sediment the flow and 
because the disturbances from the pump, which 

close the meter. 
The head lost the flume was measured four 


pressure taps connected 9.8 ft, 18.8 ft, 28.25 ft, and 
33.25 from the inlet. Settling pots were included 


Investigations Suspended Sediment Samplers,” Report No. Engi istri 
Suboffice, St. Paul, Minn., December, 1941. 
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all pressure lines trap sediment and prevent from lodging the lines. 
The static head could read the nearest 0.05 in. The temperature the 
fluid was measured thermometer the end the flume. 

was very helpful that the same sand which was used Mr. was 
still available. The sands with sedimentation diameters 0.10 and 0.16 
were used this case. Their size distributions were checked sieve 
analyses and were found Therefore, his measured values 
the settling velocity were used the present work. The size distribution 
given Fig. 


Lower Part Flow 
Upper Part Flow 


5.5 


Percentage Finer Than 
Velocity Feet per Second 


0.2 0.4 0.6 
u-Scale, per Sec 


0.088 0.124 0.175 0.246 0.01 0.02 0.03 0.05 0.10 


complete series experiments was made for each the two 
sand sizes. series had two variables, the total weight sand that existed 
the flume and the velocity the flow. Beginning with clear water, 
sand were added; the weight added each time was about the same that al- 
ready the system, following more less the terms geometric series. 
general, for each sand content, experiments were made with four flow velocities. 
They were about 1.7 per sec, 2.6 per sec, 3.5 per sec, and 4.6 per sec. 
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sample the basic experimental data given Figs. 
gives the semilogarithmic plot the velocity profiles along the center line the 
channel; Fig. gives log versus log and Fig. gives log versus The 
results obtained from all the experiments are collected and summarized 


Table 
Lower Half Flow 


100 


Number 


Values 


0.1 0.2 2.0 


noted that, except for the study the central part the channel, 
the flow was divided into two parts, one part covered nearly the upper half and 
the other part covered the lower half. part was treated alone were 
completely separate flow. This was studied and proved true 
Einstein," Assoc. ASCE, divided the section open 
trapezoidal channels into three parts, each corresponding one the sides, 
and assumed that the roughness each wall would control the flow the 
corresponding part. 


Investigation the Roughness Problem,” Schlichting, Proceedings, ASCE, 
November, 1937, following 1834. 
for the Transportation Bed Load,” Einstein, Transactions, ASCE, Vol. 107, 
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TABLE 


(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
1.78 0.00244 15.9 1.26 0.0883 
2.45 0.00452 18.0 1.21 0.12 
3.00 0.00687 17.9 1.30 0.148 
3.66 0.01026 18.0 1.28 0.181 
4.38 0.013 18.0 1.20 0.204 
5.08 0.0172 19.5 1.16 0.234 
1.78 0.00257 21.4 1.15 0.0906 
2.41 0.0045 15.6 1.26 0.12 
3.06 0.0067 21.7 1.10 0.146 
3.70 0.0093 22.0 1.10 0.172 
4.40 0.0128 20.7 1.11 0.202 
1.17 0.0011 100 21.3 1.10 0.0593 
2.01 0.00312 100 25.2 1.05 0.0998 
2.37 0.00404 0.015 21.0 1.11 0.1135 
2.39 0.00413 100 22.8 1.08 0.115 
3.53 0.00863 0.025 17.0 1.22 0.116 
1.74 0.00248 16.0 1.26 0.089 
2.60 0.00486 150 20.0 1.15 0.1246 
3.58 0.00881 0.033 20.7 1.14 0.168 
4.93 0.0155 0.057 20.8 1.14 0.222 
1.12 0.0011 19.0 1.16 0.0593 
2.40 0.0045 200 0.045 20.9 1.12 0.120 
3.29 0.00744 200 0.07 19.0 1.16 0.154 
4.23 0.0115 200 0.05 21.8 1.10 0.1915 
5.28 0.0179 0.057 20.5 1.13 0.239 
1.76 0.00239 0.02 22.5 1.07 0.0873 
2.39 0.00422 0.035 18.5 1.18 0.116 
3.26 0.0078 0.06 18.7 1.18 0.158 
1.73 0.00238 0.02 20.5 1.13 0.0872 
2.65 0.00504 0.06 17.6 1.21 0.127 
3.58 0.00835 0.08 21.6 1.10 0.163 
1.73 0.00229 0.035 21.0 1.12 0.0267 0.0855 
2.60 0.00449 490 0.08 17.5 1.22 0.025 0.12 
3.60 0.0085 0.115 18.7 1.19 0.0255 0.165 
4.89 0.0147 0.16 20.0 1.15 0.0262 0.216 
1.08 0.00101 600 18.0 0.0567 0.13 0.0554 0.057 
1.56 0.00193 0.02 18.2 1.2 0.0252 0.0785 0.124 
1.72 0.00229 0.015 20.6 1.13 0.0265 0.0855 0.124 
2.37 0.00413 600 0.06 1.12 0.0267 0.1148 0.124 
2.96 0.00582 0.12 23.4 1.05 0.028 0.1363 0.124 
3.59 0.0084 600 0.15 24.4 1.02 0.0285 0.163 0.124 
4.69 0.0142 0.17 23.7 1.05 0.0281 0.212 0.124 
1.18 0.00116 19.2 1.16 0.0608 0.122 0.0612 0.0604 
1.83 0.0025 800 19.2 1.16 0.0893 0.124 
1.75 0.00227 800 0.055 24.8 1.02 0.0287 0.0851 0.124 
2.65 0.00499 800 0.165 21.2 1.14 0.0268 0.126 0.124 
3.58 0.0086 0.26 21.0 1.12 0.0267 0.165 0.124 
4.91 0.0148 0.34 25.5 1.01 0.0291 0.217 0.124 
1.71 0.00224 0.175 25.8 1.01 0.0292 0.0845 0.124 
2.58 0.00457 0.47 20.5 1.14 0.0264 0.121 0.122 0.122 0.120 
3.59 0.00832 0.60 23.8 1.04 0.0282 0.163 0.121 0.164 0.162 
4.71 0.0136 0.60 25.5 1.01 0.0291 0.208 0.124 
1.18 0.0013 2,000 20.5 1.14 0.0644 0.124 
1.42 0.00166 2,000 21.0 1.13 0.0728 0.124 
1.71 0.00226 0.53 21.8 1.11 0.0271 0.085 0.13 0.083 0.0865 
2.61 0.00467 1.23 21.5 1.11 0.027 0.122 0.124 
3.55 0.0083 2,000 1.46 23.8 1.04 0.0282 0.163 0.124 
4.67 0.0138 1.54 21.6 0.027 0.210 0.124 


Dunes and streaks. Not able 


| 

| 
(1: 
0.35 
0.37 
0.36 
0.37 
0.37 
0.36 
0.37 
0.35 
0.37 
0.36 
0.35 
0.35 
0.35 
0.36 
0.36 
0.37 
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0.37 
0.37 
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0.35 
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0.34 
0.35 
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UNIVERSAL CONSTANT 


Top, 


(13) 


Bottom, 


(14) 


Aver- 
age, 
Kave 


(15) 


Reynolds 
number, 


(16) 


Friction 
coefficient, 
102 


(17) 


(18) 


(19) 


(20) 


(21) 


0.37 

0.356 
0.312 
0.323 
0.338 
0.348 


4.35 
6.50 


a 


ao 


a 


oun 


OFNANS 


WOW 


Bed 

(22) (23) (1) 

|Smooth 
Smooth 
Streaks 
Streaks 
Streaks 

Streaks 
|Streaks 
Streaks 


Streaks 


Streaks 
Streaks 
Dunes 

Streaks 


Streaks 


Dunes 
Streaks 


Streaks 
Streaks 
Smooth 


Streaks 

0.195 
0.197 
0.204 


see the condition the bed. 
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0.357 0.348 0.353 0.65 0.67 0.467 
0.357 0.346 0.352 0.575 0.593 0.40 0.40 
0.362 0.342 0.351 0.486 0.514 0.355 0.356 
0.366 0.346 0.356 0.366 0.388 0.198 0.262 
0.363 0.338 0.351 0.93 0.996 0.75 
0.378 0.367 0.372 0.355 0.366 0.227 0.269 
0.360 0.309 0.335 0.96 1.12 0.58 0.58 
0.375 0.320 0.348 0.577 0.69 0.438 0.438 
0.376 0.362 0.369 0.458 0.483 0.355 0.38 
0.357 0.355 0.356 0.392 0.394 0.30 0.30 
0.335 0.324 1.003 0.65 
0.347 0.335 0.638 0.685 0.437 0.51 
0.354 0.346 0.489 0.512 0.37 0.385 
0.348 0.348 0.37 0.37 0.262 0.293 
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TABLE 
FRICTION 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
1.74 0.0024 0.92 24.5 1.03 0.0285 0.0875 0.128 0.086 0.0883 
2.61 0.00476 2.90 19.0 1.18 0.0257 0.123 0.113 0.128 0.1167 
3.48 0.00818 3,600} 3.08 21.6 1.11 0.0271 0.1617 0.12 0.164 0.158 
4.89 0.0153 3.10 20.7 1.14 0.0265 0.221 0.124 
1.66 0.0025 6,800 1.49 20.5 1.14 0.0265 0.0893 0.13 0.0873 0.095 
2.61 0.00471 6,800 5.40 24.0 1.04 0.0283 0.1225 0.119 0.1248 0.117 
3.56 0.00827 6,800 5.77 23.0 1.09 0.0278 0.162 0.118 0.166 0.155 
4.74 0.0141 6.68 25.0 1.02 0.0288 0.212 0.130 0.207 0.217 
1.69 0.00277 2.39 21.0 0.0267 0.094 0.136 0.0895 0.096 
2.39 0.00415 7.23 25.0 1.02 0.0288 0.115 0.115 0.111 
2.96 0.00591 12.57 21.6 0.027 0.1373 0.112 0.144 0.126 
3.55 0.00827 12.65 24.5 1.025 0.0285 0.1622 0.105 0.174 0.145 
4.18 0.01115 12.28 26.2 1.00 0.0294 0.188 0.114 0.195 0.176 
4.91 0.01487 12.38 25.4 1.01 0.029 0.218 0.12 0.221 0.211 
1.61 0.0029 2.92 24.1 1.03 0.0284 0.0963 0.152 0.085 0.1066 
2.64 0.00535 13.66 1.07 0.0277 0.131 0.124 
3.62 0.00815 31.10 25.5 1.01 0.0291 0.161 0.116 0.166 0.156 
4.93 0.015 23,200 23.0 26.5 0.99 0.0296 0.219 0.116 0.266 0.212 
2.64 0.00471 26.3 1.00 0.1227 
1.71 0.00199 100 24.3 1.03 0.0797 
2.60 0.00476 100 1.06 0.1232 
3.55 0.00808 100 24.7 1.02 0.1605 
4.72 0.0136 100 26.5 1.00 0.2084 
1.74 0.00222 0.004 23.3 1.05 0.0842 
3.58 0.00823 0.125 24.8 1.06 0.1622 
4.60 0.01327 0.142 22.0 1.09 0.206 
1.715 0.00226 0.034 1.09 0.0584 0.085 
2.59 0.00489 700| 0.27 1.13 0.0573 0.125 
3.56 0.00823 700 0.35 24.7 1.02 0.0612 0.1622 
4.61 0.01322 0.44 23.0 1.06 0.0595 0.2058 
1.73 0.00235 1,500} 0.215 25.5 1.01 0.0622 0.0866 
100 2.60 0.00484 0.868 22.8 1.07 0.0592 0.1243 
101 3.52 0.0084 1,500 1.067 21.0 1.13 0.0572 0.164 
102 4.57 0.0131 1,500 1.05 24.6 1.02 0.0611 0.2045 
103 1.72 0.00244 0.385 25.0 1.02 0.0616 0.0884 
104 2.59 0.00494 3,100 1.94 23.2 1.05 0.0596 0.1257 
105 3.55 0.0085 3,100 2.21 21.5 1.11 0.0578 0.1648 0.121 0.167 0.1628 
106 4.58 0.0132 2.16 24.7 1.02 0.0612 0.206 
107 1.64 0.0026 6,300} 0.45 23.0 1.06 0.0595 0.091 0.126 0.0902 0.092 
108 2.58 0.00492 3.10 27.0 0.98 0.0638 0.1253 0.124 
109 3.54 0.0 851 6,300 5.0 23.5 1.04 0.060 0.1649 0.124 
110 4.59 0.01317 6,300 4.7 27.1 0.97 0.0638 0.205 0.119 0.209 0.202 
111 1.64 0.00274 0.57 27.0 0.98 0.0638 0.0935 0.134 0.09 0.097 
112 2.55 0.0054 5.1 24.5 1.02 0.0612 0.1313 0.127 0.1297 0.1332 
113 3.50 0.00886 12,700} 10.7 1.03 0.0607 0.168 0.124 
114 4.60 0.01385 13.8 26.8 0.99 0.0636 0.2104 0.121 0.2126 0.2085 
115 1.53 0.00318 25,500 1.3 22.5 1.088 0.0589 0.1008 0.148 0.091 0.110 
116 2.58 0.00563 25,500 4.8 0.95 0.0648 0.1341 0.128 0.1322 0.1362 
117 3.43 0.0096 25,500 13.5 1.00 0.0628 0.175 0.125 0.174 0.1762 
118 4.47 0.0135 25,500 24.0 0.95 0.0648 0.208 0.124 
119 1.44 0.0029 38,300 1.4 1.00 0.0628 0.096 0.154 0.107 
120 2.56 0.0059 4.5 26.7 0.98 0.0636 0.1372 0.128 0.1343 0.1402 
121 3.46 0.01005 13.0 1.03 0.0607 0.1794 0.126 0.1772 0.182 
122 4.44 0.01425 30.0 1.02 0.0616 0.2135 0.125 0.212 0.215 
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Friction 
coefficient, 
X 102 


Top Top 


(20) 


Bottom 
(21) 


Bed 


Bk ym Us dition 


(22) 


(23) 


9.75 


ded 


w 


1.63 


N 


(18) (19) 
0.97 1.03 
0.55 0.753 
0.476 0.574 
0.347 0.347 
0.786 0.93 
0.655 1.05 
0.452 0.634 
0.40 0.43 
0.87 1.07 
0.664 1.33 
0.52 0.99 
0.451 0.91 
0.436 0.653 
0.362 0.458 
0.95 0.936 
0.61 1.05 
0.538 0.843 
0.386 0.52 


1.32 
1.15 
0.76 0.89 
2.14 
1.48 
0.96 1.03 
0.82 0.87 
1.98 
1.57 
0.93 1.18 
0.82 0.91 
1.83 
1.74 
1.04 1.37 
0.85 1.07 
2.51 
1.68 
0.98 
0.855 1.33 
1.97 
1.74 
0.95 1.65 
0.88 1.56 
1.44 
0.92 1.78 
0.825 1.57 
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Reynolds 
(13) (14) (15) (16) (17) (1) 
0.341 0.313 0.327 6.55 1.98 0.562 0.65 0.224 
0.366 0.293 0.330 8.6 1.75 0.36 0.49 0.196 
0.347 0.299 0.323 12.2 1.68 0.318 0.372 0.20 
0.346 0.346 0.346 16.62 1.60 0.208 0.264 0.205 
0.386 0.30 0.343 5.67 2.26 0.482 0.593 0.207 
0.346 0.23 0.288 1.73 0.37 0.66 0.204 
0.370 0.283 0.327 0.30 0.39 0.203 
0.355 0.31 0.332 0.28 0.28 0.196 
0.343 0.260 0.301 0.56 0.62 Dunes 
0.365 0.195 0.28 0.395 0.86 0.210 ... 
0.360 0.216 0.288 0.30 0.655 0.174 ... 
0.364 0.216 0.29 0.284 0.486 0.178 ... 
0.346 0.256 0.301 0.265 0.40 0.188 
0.356 0.294 0.325 0.26 0.32 ... 
0.350 0.285 0.318 0.61 0.62 0.226 
0.326 0.221 0.274 0.263 0.56 0.181 
0.340 0.268 0.304 0.254 0.35 0.181 
0.381 0.324 0.353 0.86 
0.338 0.333 0.336 0.242 
0.349 0.322 0.336 0.250 100 
0.362 0.337 0.350 0.83 0.226 101 
0.370 0.311 0.341 0.25 104 
0.372 0.30 0.336 0.72 0.223 105 
0.362 0.327 0.345 0.64 0.209 106 
0.348 0.265 0.307 0.31 0.236 109 
0.359 0.295 0.327 0.58 0.222 110 
0.385 0.262 0.323 |Dunes 111 
0.368 0.21 0.289 0.67 0.23 
0.388 0.272 0.330 0.54 |Dunes 115 
0.380 0.215 0.297 0.78 0.208 117 
0.352 0.199 0.276 0.70 0.215 118 
0.329 0.221 0.275 1.20 1.83 0.215 120 
0.374 0.187 0.280 0.82 1.37 0.193 121 
0.354 0.183 0.268 0.60 1.2 0.20 
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Like the objective the experiments, the results can divided into two 
parts. third part will given the observations the bed. 


EFFECT SUSPENDED SEDIMENT THE CHARACTERISTICS THE FLOW 


The characteristics the flow can represented three quantities: (a) 
von universal constant (b) the momentum transfer coefficient 
and (c) the coefficient friction 


x 


Vertical Distance from the Boundary, Fee 


0.1 0.2 100 


CONCENTRATION PROFILES 


von Universal Constant, constant can calculated 


from the velocity defect law, Eq. follows: that is, 


the slope the semilogarithmic velocity profile will give 


2.3 slope the semilogarithmic velocity profile...... (9) 
but 


which equals the shear the point intersection the profile with the 
using the experimental data Bazin, Garbis showed 
that the value the shear changes along the bottom rectangular channels. 
found that the values the center the bed, divided the aver- 
age varied from 1.00 1.07 for moderate roughness. did not give any 


Turbulent Flow Open Channels,” Garbis Keulegan, Journal Research, National 
Bureau Standards, Dept. Commerce, Vol. 21, 1938, pp. 707-741. 
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data for smooth channels. few measurements, the present work, were 
made covering the entire section—runs and for clear water and run 
for loaded water. Runs and showed that was uniform along the 
boundary, whereas run the ratio the center was 
distribution interesting problem and needs more study. This dis- 
tribution might have been responsible for the small variation the values 
found for clear water channels different shapes. However, was clear 
that the distribution could contribute the variation the slopes the 
semilogarithmic velocity profiles, and, hence could introduce errors into the 
proposed method computing effect, however, seemed very small 
with respect the values observed here decreased much 50%). 
the present work was assumed uniform along 
each wall. 

find the shear velocity top and bottom and 
the flow was assumed two dimensional and, 
therefore, dr/dy which the pressure in- 
tensity and the horizontal distance along the axis 
the channel. Assuming uniform pressure gradientin 
the direction flow, will linearly related (Fig. 
depth from top bottom boundary; the intensity 
shear the top wall; and the intensity shear 
the bottom (The symbol the coordinate 7.5 
the point zero shear the point maximum velocity.) 
shear over the entire boundary. Therefore, 

For concentration per liter, the ratio the density for loaded water 
that clear water was 1.02 and its square root was 1.01; that is, would 
cause error 1%, which within the experimental error. For concentra- 
tions less than per liter the correction for the density was neglected; that is, 

These values are given Cols. and 12, Table values which 
were computed from these values and the slopes the (u- versus-log y)-lines, 
are given Cols. and both sands decreased from 0.373 about 
0.20 the average concentration increased. 

Fig. shows the relation between the values and the total load sand, 
the system for two mean velocities the channel. The points were 
somewhat scattered but Fig. shows the definite tendency decrease with 
the increase the load. the same total loads Fig. shows that the 0.16- 
sand had more effect the bottom than the 0.10-mm 
not true because, for the same total load and for the same speed, concentration 
the bottom higher for the 0.16-mm sand than for the 0.10-mm sand; and 
can assumed that both sands have nearly the same effect 

Moment Transfer Coefficient the semilogarithmic plots 
velocity profiles are given The experimental points follow reasonably 
well the logarithmic law over most the section. This proves that still 
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follows Eq. for all those different concentrations. The only changes 

The central region should excluded from this rule. Messrs. Brooks and 
Berggren’ showed that very slight change fitting the Nikuradse experiments 
will cause appreciable difference the values For example, con- 
stant will give second-degree parabola for the velocity profile, and this 
parabola can fit the velocities reported Mr. Nikuradse for 43,400 
within 0.6% over the middle half the section. evaluate 


this region from velocity measurements, because and both 


and du/dy approach zero the center. For this reason, made 
this stage evaluate from the velocity profiles. will evaluated in- 
directly studying and the relation between and 


0.40 


Total Load Sand Channel 


can observed from Fig. that there was deviation from the logarith- 
mic law high concentrations near the bed the stream. This occurred 
Mr. work. 

Coefficient Friction, furnished the relation between the 
friction coefficient and the Reynolds number for smooth rectangular 
channel: 


value the parameter y/v, which the transition between the boundary 
laminar film and the inner turbulent region corresponds. 
Taking 11.6 according the experiments Mr. Nikuradse, and 


Fig. shows such curve with the experimental results. The curve fitted 
the points the clear water quite satisfactorily. This proves that the walls 


0.30 Region Sand Size 
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LEGEND 
Symbol Run 


0.025 
0.020 


~ 


0.014 
0.040 


LEGEND 


Run 
0.035 


Values the Friction Coefficient, (Multiply 107) 


0.017 


0.16 SAND 


0.014 


Load 
Kilograms 


23.2 


Load 
Kilograms 


122 119 38.3 

115 118 25.5 

111 116 12.7 

0.030 107 110 
103 106 3.1 

102 


0.030 
13.2 
3.6 
<2.0 
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are smooth. The general pattern the results plotted Fig. agrees very 
well with the results collected Mr. O’Brien and 

Taking the curve versus Eq. reference, can observed that, 
for each concentration each sand, there was certain critical Reynolds 
number below which the friction coefficient began increasing with respect 
that the clear water, and above which the coefficient friction for all con- 
centrations was practically the same. was observed from the glass windows 
that the critical value corresponded the velocities that were able 


2.0 
LEGEND 
Upper Portion 
0.4 Lower Portion 0.8 
0.2 0.4 0.6 0.8 1.0 0.4 0.8 1.2 1.6 2.0 


sweep all the dunes from the floor. Below these critical values appre- 
ciable fraction the sand was moving bed load dune forms, and above 
these values there were dunes and most the sand was carried 
suspension clouds. This fact was also observed from the values given 
Col. Table and from Above these critical values was difficult 
increase since most the sand was already suspension. just re- 
distributed itself, shown the crossing and overlapping the concentration 
profiles (Fig. 6). 


SEDIMENT DISTRIBUTION AND SEDIMENT TRANSFER COEFFICIENT 
The two sets curves given Figs. and show the sediment distribution 
over the section. Fig. was drawn verify Eqs. whereas Fig. was drawn 


verify Eq. 
For most the depth the channel (excluding the central part) Fig. 
shows that Eqs. hold well. The slight curvature the top wall will 


Engineering, Univ. California, Berkeley, Calif., 1937, Vol. No. pp. 343-383. 


2.4 

oO 
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discussed subsequently prove that this deviation due sand grading. 
According the slope the lines log and log must equal 
The values measured this method are given Cols. and 21, Table 
whereas the calculated values w/k are given Cols. and 19. The 
upper and lower halves the channel were treated separately. Fig. shows 
the relation between the measured values (which should 
and the values w/(k The linear relation both curves for the two sands 
shows that the assumption Eq. fairly reasonable. For the 0.10-mm sand, 
1.5; and, for the 0.16-mm sand, 1.3. 

Fig. helps studying the central part the flume. The straight parts 
the curves (which always occur the center) indicate that the assumption 
constant here quite satisfactory. The constant value the central 
parts the flow (which will called can evaluated directly from Eq. 
and from these curves. The nondimensional values are given 
Col. 22, This value nearly constant for most the runs, for both 
the sands, and its average value about only runs that are far from 


Error Curve 
and 


this value are those low speed and high concentration. the latter runs big 
dunes were observed and the conditions the flow were believed unsteady 
because the periodical disturbance caused the travel dunes past the 
measuring point, these readings were excluded. The straight part, and 
accordingly the constant was observed cover least the middle third 
the channel depth each run (Fig. was measured and compared 
with the parabolic form Eq. was found that would intersect the 
third points, asin According the present experiment, Fig. would 
the most reasonable approximation the dimensionless form 
previously mentioned, could also represent the form €m. 


Loap 


The bed load was observed through the glass windows. For high speeds and 
high concentrations was not possible see what was occurring. For small 
loads sand and for low velocities the bed load was the form streaks. 
These streaks were not very steady and their particles were moving steps. 
When the velocity was increased, all the streaks disappeared and the particles 


went into suspension. 
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When the concentrations were somewhat increased for low speed, dunes 
began appear the bed isolated individual dunes crescent shape, small 
size compared the channel width. the beginning, single dune might 
appear every along the bed. When the concentration was increased, 
the number those dunes increased and they started joining each other form 
one dune that extended across the channel. The form the bed load de- 
scribed letters Col. 23, Table 

was observed for the maximum loads reached the 
case the 0.16-mm sand. During run 119, which was made with the low 
velocity 1.7 per sec, very large dunes waves were observed the bed. 
They occurred irregularly every now and then; sometimes two might follow 
each other closely and sometimes more than half hour elapsed between these 
waves. wave was about in. long and 0.6 in. maximum height and 
would cover the full width the flume. The dunes and waves sand that 
run were large and unsteady that the entire flow was considered unsteady; 
then the results the run were rejected. 

With the same concentration but with higher velocities seemed the 
dunes and waves were smoothed out and that uniform layer was moving along 
the floor with small dunes its surface. run 120 this layer was about 0.25 
in. thick, and the dunes were about 0.05 in. high over the layer. 


RESULTS 
Effect Suspended Load the results the studies can 
appraised discussion the effect the suspended load four character- 
istics the flow: (a) the universal constant; (b) the momentum transfer 
cient; (c) the friction factor; and (d) the velocity distribution. 


(a) The Universal Constant, his similarity theory turbulence 
pattern Mr. von introduced the factor proportionality be- 


du/dy 


tween the velocity derivatives and the mixing length 


the Prandtl mixing length. 

Assuming that the turbulence pattern will remain similar from point 
point, will remain constant over the cross section the flow. The decrease 
the value which was observed Mr. Vanoni, varied from run run, 
but each run was the same for the entire profile. present work, for 
the upper part the section was constant; and for the lower part was also 
constant but had different value than for the upper part. This means that the 
change was not local from point point following the different concentra- 
tions these points, but changed from section section and remained constant 
over each section. seems, then, that the presence certain proportion 
sediment suspension will affect the main pattern the turbulence the 
entire region from the controlling boundary the zero shear plane. 

The slight curvature the semilogarithmic velocity profiles observed the 
case high concentration near the bed (Fig. cannot explained with any 
local affected the local concentration, should decrease 
the bed approached and this will curve the line the opposite direction 
that shown Fig. This curvature was observed occur only near the 
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bottom for appreciable concentrations and increased with the increase the 
load. This might evidence that the curvature due some interaction 
between the flow and the presence the particles the bed. The question 
the local effect sand has been raised Van 

measuring the velocity profiles the wind above circular concrete 
pavement, observed approximate variation the values 
from 0.6 the lower atmosphere. The constant turbulent transfer 
decreases the atmospheric stability increases. Mr. Sheppard declared these 
values are very rough; the wind profiles the experimental site were not repre- 
sentative the boundary layer because the limited extent the measuring 
surface (concrete.) 

The only source error evaluating can due large variation the 
along the boundary. Run showed that the existence sand 
could never change the take care even small part the 
changes the semilogarithmic velocity profile slopes. 

Momentum Transfer Coefficient, has been shown how was 
and how Fig. represents its distribution. The relation 
between and was verified the part the section where both and 
could measured directly. very reasonable assume that this relation 
will hold over the remainder the section—that is, the central part, where 
isnot known. Therefore, Fig. can represent the distribution well 
This curve agrees very well with the curve given Messrs. Brooks and 
Berggren’ under their assumption the error curve (see Fig. present 
relation easier handle. 

(c) The Friction Factor, results the present work given 
Fig. agree very well with the results reported Messrs. O’Brien and 
graphically (as Josef agree with the other 
results. Below certain Reynolds numbers increases with the concentration, 
whereas, for greater values the concentration has effect. Comparing the 
different curves versus given this paper and clear 
that the curves the different concentrations approach the clear water curve 
different ways which depend mainly the size the suspended material. 
Messrs. O’Brien and show how the resistance the clay slurry 
deviates sharply from the clear-water curve. This similar Fig. 10(a) for 
the 0.10-mm sand. When the size the sand was increased 0.16 (Fig. 
the curves approach the clear water more smoothly over longer region. 

From the careful observation the bed the present work, these critical 
velocities are related the existence dunes the clear that 
for the same load and velocity the bed load will increase with the size the 


Vanoni, Transactions, ASCE, Vol. 111, 1946, 106. 

Aerodynamic Drag the Earth’s Surface and the Value von Constant the 
Lower Sheppard, Proceedings, Royal Soc. London, January, 1947, 208. 

Vol. 104, 1939, pp. 1334-1380. 

Discussion Josef Montgomery ‘‘Transportation Sand and Gravel Four-Inch 
Howard, ASCE, Vol. 104, 1939, 1350, Fig. 

Engineering, Univ. California, Berkeley, Calif., 1937, Vol. No. 361, Fig. 
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suspended material. also clear that the chances forming dunes increase 
with the increase the bed load. This the reason for the dependence the 
critical Reynolds numbers the size the suspended material. The initial 
roughness might affect these critical numbers slightly. channel has very 
rough walls, say, the case Fig. 10, can expected that the reference line 
clear water will raised. There will still transition region, but the 
critical values might reached smaller values than the case 
smooth wall. 
With the indicated substitutions, Eq. for smooth channels becomes 


which represented Reynolds-number form for the thickness the 
laminar sublayer follows: was found empirically that 
11.6 for the clear fluids. For rough walls follows the relation: 


which the equivalent roughness height parameter and the wall 
velocity. 

Two main factors will affect factor the internal system flow 
represented and can called the internal factor. The other one the 
boundary conditions represented either and the walls are smooth (or 
the roughness height they are rough); and this can called the boundary 
factor. There are some other factors, but they appear terms secondary 
order magnitude and can neglected. 

Varying the internal factor varying alone and holding the boundary 
factor constant, will change approximately the square according 
Eq. 14. must decrease the presence the suspended material, because 
the decrease for the effect the boundary factor, can seen that 
with sediment-laden flow and continuous interchange particles between 
the bed and the fluid, quite difficult visualize the picture the boundary. 
conceivable that the laminar sublayer will broken the continuous 
bombardment the sand particles through it. This will increase and the 
bottom region will somewhat similar the case near rough surface even 
looks smooth. When there are dunes the bed, obvious that they 
greatly increase the bed roughness and will increase. Now when both the 
factors are considered, will decrease increase according the ratio between 
the effects the two factors. 

low velocities, below the critical value, the dunes will increase the bound- 
ary effects and accordingly will increase, seen from Just after the 
dunes are swept off the bed, the case the fine particles, seems the 
effect exceeds the boundary-layer effect the sediment the bottom and 
less than for the clear water. This can seen Fig. 10(a). high 
velocities the two effects are nearly balanced and not changed. 
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(d) Velocity shows some the velocity profiles 
the center line, for high concentrations, compared the profiles clear water 


flow, which have nearly 
the same values 

Fig. 13(a) shows the 
effect dunes the 
velocity distribution. 
shows how the maximum 
velocity and, correspond- 
ingly, the plane zero 
shear were shifted upward. 
This indicated that the 
shear the bed was in- 
creased (Fig. causing 
increase with the dunes. 
When the effect the 
dunes was not strong 
the other factors came into 
the picture and resulted 
profiles such Figs. 
13(c), and 13(d). 

For high velocity was 
shown Fig. that 
(or u/U;) was the same for 
both clear and sediment- 
laden water. Fig. 13(d) 
and Col. 10, Table 
showed that these cases 
the maximum velocity al- 
ways lies the center; 
therefore, the top 
must equal the 
bottom. There 
other way explain the 
big difference between the 
velocity profile slope 
the top halfand the bottom 
half except the change 
compare the two sides 
the velocity profile for 
high concentration (see 
Fig. 13(d)). velocity 
near the bottom shows 
quite decrease the wall 


U=1.53 


Run U=2.45 
Run 116; U=2.58 


0.6 
1.2 


Values 


U=3.66 
U=3.43 


0.05 0.10 0.10 
Feet 


0.05 


the roughness that wall had been increased. This could explained the 
effect the sand breaking the laminar sublayer and decreasing its 
thickness the same way increasing its roughness. 
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Sediment Distribution—The sediment concentration curves Fig. are 


slightly curved near the upper surface the channel. The concentrations 
this region were usually very low; and, when some the samples were examined 
carefully eye, they appeared finer than the sand customarily used. 
known that however carefully the sand sieved will always have certain 
standard deviation particle sizes, shown Fig. the sample 


Um 


Vertical Depth from the Boundary, 


Curve Based 
Shown Fig. 


0.2 


0.4 


0.6 


Center Line 


Cc 


THE 


into different sizes, clear that the very fine portion will more uniformly 
distributed than the coarser portions. Now, knowing that Fig. for log and 
that the concentration very low the top, clear that the small proportion 
the very fine uniformly distibuted particles will amount appreciable 
fraction the sample the top, whereas negligible the bottom. This 
grading then the cause the observed curvature. 


ahs 

0.6 

0.4 
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was shown the results that 1.5 for one sand and 1.3 
for the other sand. These two observed values are not sufficient establish 
the relation between and the size sand. More sizes and mixtures sizes 
should tested the future. Messrs. Sherwood and Woertz‘ found the eddy 
diffusion the main body the turbulent fluid for the case vaporization 
water into turbulent gas stream equal 1.6 €m. found 


LEGEND 
Runs 61, 62, 63, 
Runs 65, 66, 67, 
Runs 79, 80, 81, 
Runs 101,-102 
Runs 105, 106 
Runs 112, 113, 114 
Runs 116, 117, 118 


Curve Based 
Shown Fig. 


Vertical Depth from the Boundary, 


SEDIMENT CONCENTRATION ym) 


that the heat transfer coefficient between two parallel plates, varied between 
1.5 and 1.22 with average 1.3 

Fig. represents the assumed distribution the nondimensional form 
across the central profile channel. Fig. shows curve drawn 


Gradients Turbulent Air Corcoran, thesis presented the Califor- 
nia Inst. Technology Pasadena, 1948, partial fulfilment the requirements for the degree 
Doctor Philosophy. 
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for the sediment concentration distribution under this assumption. must 
noticed that the measured values were always used for plotting the points 
this curve. The points twelve runs for the 0.10-mm sand and ten runs for 
the 0.16-mm sand are shown Fig. 14. The agreement between the experi- 
mental results,and the assumed theoretical curve remarkable, which proves 
that the assumption Fig. reasonably correct. 


Top 
1.353; 34 0.94 0.94 0.90 0.92 1 
1.24 1.55 
148¢e o15 024 
2.0 
230 2.2 
480 7 45 04.0 04.0 045 
v 
8.0 5 = 
21.7 
0245 27028 27.7 021.0 26.8 “35 
38.1 47.0 52.3 49.3 53.0 
Bottom 


CONCENTRATION DISTRIBUTION OVER THE SECTION FOR No. 
(GRAMS PER LITER 


Fig. shows the concentration distribution over the entire cross section 
71. clear that the middle half the contours run parallel the 
top and bottom, showing that the assumption two-dimensional flow very 
reasonable. Near the side walls the concentration more uniform over 
vertical section, which indicates that stronger near the side walls than 
the center. Mr. Van has explained the more uniform sediment con- 
centration near the wall saying that any vertical line more uniform 
near the wall. This could not the reason because, seen from Fig. the 
average slope the concentration profile was mainly controlled the average 
value the profile and not its variation. 


The von constant decreases when the load suspension 
increased. Both the sands used the experiments have nearly the same 
effect which decreased low 0.20 for load sand equal kg., 
other words, when the sand content the system was per liter. The 
decrease indicates damping the turbulence. 

The change does not follow the change the concentration from 
point point over one cross section, but varies from section section, 
maintaining constant value over each one. This value decreases with the 
increase the total quantity suspended sand the section. 

The value the momentum transfer coefficient affected the 
presence the sand only through the changesink. The form its distribution 
curve will remain the same for sediment-laden water for clear water. This 
coefficient cannot measured directly the center the channel but, because 
similar €,, one can reasonably assume that follows Fig. 12. 
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The coefficient friction for stream carrying suspended sediment 
exceeds that for clear flow only when dunes form the bed. there are 
dunes, the present experiments show that can hardly affected the 
presence sand. This for smooth walls and probably for walls roughened 
particles smaller than, equal to, the suspended particles. 

The sediment transfer coefficient follows approximately 
equal 1.5 for the 0.10-mm sand and 1.3 for the 0.16-mm.sand. 

Two devices—the basic equation for sediment distribution Eq. and the 
distribution curve vertical profile the channel (see Fig. 12)—can 
define, accurately, the concentration distribution represented Fig. 14. 
The only difficulty this procedure that and cannot predicted from the 
theory and they must assumed determined experimentally. 
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APPENDIX. NOTATION 


The following letter symbols, introduced this paper, are assembled for the 
convenience discussers: 
concentration suspended material, expressed mass per unit 
volume; denotes value arbitrary reference level; 
concentration over given cross section; denotes measured mean 
value 
diameter sand particles: 
geometric mean sieve size; 
sedimentation diameter sand; 
equivalent roughness-height parameter; 
acceleration gravity; 


_ = 2-1; 


hydraulic gradient slope; 

von universal constant turbulent exchange: 
the constant the bottom region section; 
the constant the top region section; 

mixing length; 


the characteristic value the parameter 
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pressure intensity; 

Reynolds number; 

hydraulic radius; 

temperature; 

velocity flow; mean velocity over the cross section: 


the Prandtl friction velocity 


value the top wall boundary; 


velocity (see also U): 
turbulent velocity; 


average velocity along the center-line profile; 


Umax Maximum value 
total load sand the flume channel; 
settling velocity particle still fluid; 
horizontal distance along the axis the channel; 
vertical distance from the boundary: 
measured from the boundary the plane maximum 


velocity, also, the coordinate the point zero shear; 


the distance from the top boundary the bottom bound- 


factor proportionality; 

thickness laminar sublayer; 
transfer coefficient: 


ary; 


momentum transfer coefficient; 
sediment transfer coefficient; 


friction coefficient; 

kinematic viscosity; 

density fluid; mass per unit volume; 

standard deviation; geometric standard deviation; and 
shearing stress: 


value the point intersection the profile with the 
boundary; 

value the bottom boundary; 

value the top boundary. 
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